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.  A  PROGRAM  OF  AUTOMATIC  DIFFEREMTATION 
FOR  THE  BESM  COMPUTER 

'  [This  is  a  translation  of  a  booklet  by  L.  M.  Vega, 

L.  N.  Korolev,  N.  V,  Sukhikh,  and  T.  S.  Frolova, 

Institute  of  Precision  Mechanics  and  Calculating 
Techniques,  Moscow,  1939*3 

The  paper  contains  a  general  description  of  a  pre¬ 
gram  and  a  method  of  analytical  differentiation  for 
high-speed  electronic  digital  computer  BESM. 

Some  results  of  this  program  application  to  con¬ 
crete  mathematical  expressions  are  presented  and 
the  program  block-diagrams  are  given. 

The  obtaining  of  derivatives  of  complex  mathematical  expres¬ 
sions  is  a  rather  laborious  process.  The  task  of  obtaining  the 
derivatives  by  means  of  a  computer  is  a  problem  interesting  from 
the  point  of  view  of  programming  logical  problems;  it'  also  has  an 
applied  significance. 

In  this  work  there  is  expounded  one  of  the  methods  of  auto¬ 
matic  differentiation  of  mathematical  expressions  and  there  is  des¬ 
cribed  the  arrangement  of  a  corresponding  program,  developed  by  a 
group  of  employees  of  the  Institute  of  Precision  Mechanics  and 
Calculating  Techniques,  Academy  of  Sciences  USSR  for  the  BESM 
computer. 

Let  us  agree  that  by  an  expression  we  designate  a  line 
consisting  of  letters,  numbers,  symbols  of  arithmetical  operations, 
symbols  of  elementary  functions  and  the  symbols  of  open  and  closed 
brackets.  In  order  for  an  expression  to  have  mathematical  signifi¬ 
cance,  the  symbols  in  it  must  combine  in  a  definite  way;  for 
example,  the  number  of  open  brackets  must  always  be  equal  to  the 
number  of  closed,  also  two  symbols  of  multiplication  one  after  the 
other  should  not  occur,  etc.  Therefore,  naturally,  we  will  deal 
only  with  expressions  which  have  mathematical  significance. 

Let  us  call  parameters  or  arguments  (or  simply  variables) 
the  letters  and  numbers,  i.e.,  all  symbols  except  brackets,  symbols 
of  arithmetical  operations,  symbols  of  elementary  functions.  Any 
expression  may  be  regarded  as  a  function  of  the  parameters  of 
which  it  is  composed.  One  and  the  same  function  may  be  given  in 
various  ways.  For  example,  the  expression  a  +  b  +  sinVx,  which 
is  obviously  a  simple  function  of  a,  b,  x  transforms  into  a  complex 
,  function  of  arguments  a,  b,  x: 
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by  means  of  the  following  set  of  equalities: 

“»-«♦»;  «a  -  «ia  «4  ;  u4  -VT. 

of  var^!ftf'er'  fY  f^33*011  can,  by  means  of  substitution 
i1  L1*  put  int0  the  trivial  form  f (u)  .  «  (where  u  is  a 

•Sod  of'difS  °r^Sinal  &r&mentB).  The  essence  of  Ly 

Mjr»«jtsUtxo«»  is  m  the  utilization  of  the  rule  fo7 

S  feJ^iati0C'  of  a  C03Kplex  function.  The  simplicity  or  complexity 

JESS'S  Sf?4 “  -  “• «  SwSSS* 

is  presented  in  the  form  of  a  complex  function. 

,  tbe  oS'Sfs  °f  the  offered  method  lies  an  algorithm  of 

introduGtif.«XI>^nSi0n  °f  ths  initial  expression  by  means  of  the 

terafarftL  W  var*ables  in  the  set  of  equalities,  whose  left 
w  M  ®  the  new  variables,  and  whose  right  are  easily  differen- 

viSiblestn^  eXpression^  iSo  which  oid  aS  n^  ■ 

Session  £  V^able  sutstltutes  for  the  entire 

E  equalities  and  equalities  obtained  by  means  of 

Sitiaf  e^rasS^  Uni*uel3r  the  derivative  of  the 

the  last  differentiatedCequali^nofl5theUnewSvariable«^iiUei0nbyn^0 
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For  example,  the  egressions  cited  below  are  pairs: 

8laa  f  Btb--  «-*  ;  «•«  etc. 

pairs  f°r  e:Cpfnsion  of  the  initial  expression  into 

lation  ^+Z  +  ?e  e^lUivalent  to  the  algorithm  of  formula  trans- 
^  tSee  SdSS1  operations  in  the  formula)  into  the  sequence 
command  in  a  thLeTS^  ^  aatomatic  programming.  The  arithmetical 

m  corresponds  to  an  «l«wrta*y  substi- 

element^^0pa^r^a^2^df??°^a^Ce^a°^^i^ei®x^^®^a^°:nt°nt°CCeSSiVe 

sssrsLsas  «. 

new  parameters  into  the  derivative  of  the  last  pair. 

separatfmS/^S10  th®  deSCrlption  of  the  arrangements  of 
P  rate  parts  of  the  program  of  automatic  differentiation,  let  us 
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analyze  ■toe  accepted  system  of  coding  of  all  toe  symbols  which  are 
found  in  it.  The  coding  system  was  conditioned  by  toe  arrangement 
of  the  BESM  computer. 

For  each  symbol'  eight  binary  columns  are  assigned,  i.e., 
in  each  cell,  beginning  with  the  32nd  column,  4  symbols  may  be 
placed. 

The  elementary  pairs  consist  either  of  three  symbols  or  of 
two;  therefore,  each  pair  takes  up  no  more  than  one  cell. 

The  parameters,  symbols  of  elementary  functions,  symbols  of 
'arithmetical  operations,  and  brackets,  are  coded  by  numbers  from 
O.i  to  5.15  (sde  Table  l). 

For  convenience  in  coding  toe  elementary  pairs  Uq,  up...^, 

to  which  the  elementary  expression  is  expanded,  those  are  taken 
which  coincide  with  the  address  of  block  cells  in  which  toe  corres¬ 
ponding  pairs  are  stored.  For  storing  elementary  pairs,  a  special 
location  of  cells  in  the  IMD  (internal  memory  device)  of  toe  com¬ 
puter,  from  6.0  to  10.15  is  assigned?  hence  the  pair  uQ  is  coded  as 

6.0,  the  pair  u^  as  6.1,  etc.  The  code  of  a  derivative  pair  or  of  a 
parameter  is  obtained  by  adding  the  code  of  the  pair  or  parameter  to 
toe  constant  5.0?  for  example,  Uq  is  coded  as  11.0,  as  H.l,  etc. 

For  storing  derivatives  of  elementary  pairs,  another  location' 
in  the  IMD  is  assigned,  from  11.0  to  25.15.  codes  of  derivative 
pairs  do  not  coincide  with  the  numbers  of  the  cells  in  which  the 
appropriate  derivatives  are  stored,  since  toe  pair  derivative  may 
take  up  several  cells. 

In  order  to  avoid  coincidence  of  toe  codes  of  parameter 
derivatives  with  pair  codes,  it  is  possible  to  construct  another 
system  of  coding,  but  that  would  complicate  the  perforation  of  for¬ 
mulas,  since  more  than  8  columns  would  be  needed  for  coding  each 
symbol. 

Each  expression  subject  to  differentiation  on  toe  computer  is 
brought  to  a  certain  standard  form.  At  insertion  into  toe  computer, 
the  expression  is  enclosed  in  brackets.  The  numbers  which  are  part 
of  it,  with  the  exception  of  0.5  and  whole  numbers  from  1  to  9 
inclusive,  are  replaced  by  letter  symbols.  The  operation  extraction 
of  a  root  is  replaced  by  raising  to  toe  power  of  0.5,  the  complex 
arguments  being  enclosed  into  brackets.  The  expression,  coded  with 
consideration  of  these  rules,  is  introduced  into  the  IMD  of  toe 
computer.  It  should  be  noted  that  toe  rules  of  coding  do  not  place 
any  limitations  on  toe  structure  of  toe  expression  and  are  induced 
by  considerations  of  programming  convenience. 

The  following  division  of  toe  computer  memory  has  been 
accepted. 

The  expression  being  differentiated  is  placed  in  toe 
location  of  cells  from  1.0  to  5*15* 

The  elementary  pairs  are  formed  in  the  location  of  cells  from 
6.0  to  10.15.  The  derivatives  of  elementary  pairs  are  obtained  in  a 
location  of  cells  from  11.0  to  25.15* 
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The  derivative  of  the  whole  expression  is  formed  in  the 
location  of  cells  from  26.0  to  30.15. 

Let  us  examine  the  first  part  of  the  program  of  expansion 
into  pairs. 

In  the  formation  of  pairs,  i.e.,  in  representing  a  mathema¬ 
tical  expression  by  means  of  a  succession  of  pairs,  the  rule  of 
arithmetical  operation  must  be  taken  into  consideration.  Therefore, 
at  first,  the  bases  of  powers  with  exponents  are  joined  into  pairs, 
the  symbols  of  elementary  functions  with  parameters;  then  pairs  of 
parameters  are  formed,  connected  with  s  ymbols  of  the  mathematical 
operations  of  multiplication  or  division;  and  then  pairs  which 
contain  the  symbols  of  addition  or  subtraction.  In  reading  a 
mathematical  expression  successively  symbol  after  symbol  from  left 
to  right,  the  symbols  +,  .,  :  serve  as  signals  of  pair  formations 

which  contain  a  power  of  a  symbol  of  an  elementary  function;  the 
symbols  +,  -,  )  serve  as  signals  for  pair  formations  Which  contain 
the  signs  .,  :;  and  ()  serves  as  a  signal  for  +,  -,  pair  formations. 

Example.  Let  the  expression  [In  sin  x  +  e  -va*x/*  arctg  x| 
be  differentiated  with  respect  to  x.  This  expression  is  expanded 
into  the  following  pairs; 

36 


6.0 

sin  x 

6.1 

in  6.0"' 

6.2 

a*x 

6.3 

0-6.2 

6.4 

eX  6.3 

6.5 

arctg  x 

6.6 

6.4*  6.5 

6.7 

6.1+6. 6 

The  signal  for  pair  formation  is  the 
The  signal  for  pair  formation  is  the 
The  signal  for  pair  formation  is  the 
The  signal  for  pair  formation  is  the 
and  the  proceeding  sign 
The  signal  for  pair  formation  is  the 
The  signal  for  pair  formation  is  the 
The  signal  for  pair  formation  is  the 
The  signal,  for  pair  formation  is  the 


sign  + 
sign  + 
sign  ) 
sign* 

sign  • 
sign  ) 
sign  } 
sign  ) 


The  process  of  representing  the  expression  in  the  form  of  a 
sequence  of  pairs  is  best  followed  by  using  the  instruction  diagram 
of  this  part  of  the  program  (see  Instruction  Diagram  l). 

For  a  negative  number  or  -ui,  if  it  is  the  indicator  of  a 
power  or  if  it  follows  directly  after  an  open  bracket,  the  corres¬ 
ponding  pair  o  -x,  o-  ui,  i.e,  o  minus  this  number  is  put. 

In  the  program  for  obtaining  a  sequence  of  pairs  from  a 
mathematical  expression,  there  are  two  instructions  for  formation 
of  pairs.  One  instruction  forms  pairs  from  two  symbols  of  the  type 
Inx,  the  other  forms  pairs  from  three  symbols,  for  example,  a*x, 
a+b,  etc.  Thus,  each  pair  is  formed  in  one  of  these  instructions. 
The  instruction  for  recording  the  pair  in  location  6.0  -  10.15  is 
a  common  one.  Consequently  each  pair  passes  through  it. 

The  instruction  for  storing  symbols  in  location  26.0-39*15 
is  also  common  for  all  pairs.  In  this  location  the  symbols  are 
recorded  one  to  each  cell,  in  a  sequence  from  the  initial  expres¬ 
sion,  which,  during  the  work,  is  "read"  once  from  left  to  right. 
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Recording  In  location  26.0-39*15  takes  place  until  the  symbol 
selected  from  .  the  initial  expression  turns  out  to  be  the  signal 
that  a  pair  may  be  obtained.  This  pair  is  placed  in  an  appropriate 
place  for  the  location  of  pairs,  and  in  location  26.0-39*15  there 
occurs  substitution  for  two  or  three  symbols  of  the  formed  pair 
by  one  symbol  --  the  number  of  the  cell  in  which  the  pair  has  been 
recorded.  later,  the  computer  either  uses  this  symbol  for  the 
formation  of  the  next  pair,  or,  depending  on  the  instruction, 
continues  to  select  the  next  symbols  of  the  original  expression, 

*  The  intermediate  recording  of  symbols  in  location  26.0-39*15 

simplifies  to  a  considerable  degree  the  program  of  algorithmic 
expansion  into  pairs,  since  it  liberates  it  from  the  complicated 
work  of  substituting  variables  within  the  initial  expression  and 
allows  it  to  move  along  in  only  one  direction. 

In  the  program  there  is  an  instruction  for  the  "destruction 
of  brackets,  which  operates  when  there  turns  out  to  be  only  one 
symbol  within  the  brackets.  If,  for  example,  we  come  across  the 
combination  (  u  +  a  ),  then,  after  formation  of  the  pair  u  +  a  and 
the  substitution  for  it  of  u^,  we  obtain  (uj_).  These  brackets  are 
de stroyed. 

A  derivative  pair  which  has  passed  through  one  of  the  above- 
mentioned  instructions  is  recorded  in  the  location  of  cells  from 
59.0  to  59.4.  The  instruction  for  storing  the  derivative  of  the 
given  pair  in  this  location  is  common  for  all  pairs.  Simplification 
of  the  obtained  expression  of  the  derivative  is  then  performed, 
since  the  differentiation  is  performed  according  to  formulas  of  a 
general  type;  for  example,  for  differentiating  xa  or  ax,  the  base 
and  the  exponent  are  substituted  into  the  general  formula  for  the 

derivative  of  a  power:  ‘ 

(uv)'  -  u^lau  ♦ 


Die  simplification  is  further  reduced  by  destroying  brackets, 
unnecessary  for  the  expression  of  the  derivative  of  a  concrete 
pair,  symbols  of  elementary  operations,  zero-addends,  unit- 
multipliers  and  products  into  which  zero  enters  as  a  factor.  For 
each  instruction  for  derivative  derivation,  there  is  a  particular 
instruction  for  simplification  of  the  derivative  (see  Instruction 
Diagram  2).  Such  simplifications  as  collecting  terms  and  reduction 
of  fractions  are  not  performed. 

The  derivatives  which  have  passed  through  the  .simplification 
instructions  are  recorded  in  the  location  of  cells  from  11.0  to 
25.15;  furthermore,  this  record  is  realized  in  general  for  all 
instruction  pairs.  If  the  derivative  of  a  pair  turns  out  to  be 
zero,  then,  in  the  corresponding  cell  of  location  11.0-25.15,  the 
code  for  zero  is  recorded.  In  the  same  instruction,  the  relative 
address  of  the  origination  of  the  derivative  according  to  11.0- 
25.15  (obviously  Nrel*=Ncell  -11.0)  is  obtained.  Biis  address  is 

recorded  in  the  last  8  columns  of  a  cell  of  location  6. 0-10. 15, 
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where  the  corresponding  elementary  pair  of  this  derivative  is 
located.  Thus,  each  elementary  pair  has  the  address  of  its 
derivative.  For  example,  pair  ug  has  the  form  a«x  in  cell  6.2 

(see  example,  page  6),  has  o«x-w.»l  in  cell  59*0  and  a  in  cell  11.3* 

As  a  result  of  treating  each  elementary  pair  by  the  second 
part  of  the  program,  a  series  of  elementary  derivatives  is 
obtained. 

As  an  .example  we  will  cite  the  contents  of  the  location  of 
^  cells  from  6.0  to  1.0.15  and  the  location  of  cells  from  11.0  to 
*25.15  which  will  be  obtained  as  a  result  of  the  work  of  the  second 
part  of  the  program  for  automatic  differentiation  (see  Table  2). 

By  Ihe  beginning  of  the  third  part  of  the  program  pairs  are 
recorded  in  the  location  6.0-10.15;  furthermore,  in  the  first 
eight  columns  of  cells  of  this  location  there  are  decreased  by 
11.6,  the  numbers  of  the  cells  in  which  the  corresponding  deriva¬ 
tives  originate,  and  in  location  11. 0-25*15  the  derivatives  of  the 
pairs  are  recorded. 

The  work  of  the  third  part  of  the  program  terminates  with 
obtaining  the  derivative  of  the  initial  mathematical  expression, 
recorded  in  the  location  26.O-30.i5. 

We  know  that  the  last  pair  replaces  the  initial  expression; 
consequently,  the  derivative  of  the  last  pair  equals  the 
derivative  of  the  initial  expression.  Therefore  it  is  natural 
to  obtain  the  derivative  of  the  initial  formula,  expressed  in 
the  initial  parameters,  by  subsequent  substitution  into  the 
derivative  of  the  last  pair  of  values  for  all  new  parameters 
which  are  a  part  of  it  and  which  signify  pairs  and  derivatives  of 
pairs. 

Let  us  describe  in  more  detail  the  process  of  obtaining 
the  derivative.  In  the  last  finished  cell  of  the  location  6.0-10. 15 
there  is,  decreased  by  11.0  (relative),  the  number  of  the  cell  in 
which  the  derivative  of  the  last  pair  originates.  The  obtained 
number  allows  selecting,  from  location  11.0-25.15,  the  derivative 
of  the  last  pair  and  bringing  it  into  location  1.0-5 .10.  From 
location  1.0-5.10  consecu.  tively  8  columns  are  taken,  which  are 
checked  for  coincidence  with  the  codes  of  the  parameters  that 
denote  pairs  or  derivatives  of  pairs.  If  it  turns  out  that  there 
is  no  such  coincidence,  the  computer  registers  the  code  of  the 
symbol  registered  in  these  eight  columns  in  location  26.0-30. 15. 
However,  if  in  the  selected  8  columns  there  is  a  parameter  code 
which  denotes  a  pair  (or  derivative  .pair),  then,  from  location 
6.0-10.15,  the  pair  corresponding  to  that  parameter  (or  from  loca¬ 
tion  11,0-25.5  the  derivative  of  the  pair)  is  selected  and  is 
transmitted  into  location  63.ll-63.i5,  where  it  is  analyzed  and 
then  recorded  in  location  26.O-3O.i5  in  brackets  or  without  brackets 
depending  on  the  results  of  the  analysis  (see  Instruction  Diagram 

5). 

Having  exhausted  all  elements  of  location  I.O-5.I5,  the 
computer  replaces  them  by  elements  of  location  26.0-30. 15  and 
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repeats  the  same  cycle*  This  continues  until.  In  location  26.0-30.15#  * 
the  derivative  of  the  initial  mathematical  expression  is  obtained, 
i.e.,  until  the  expression  recorded  in  location  26.0-30. 15  finally 
contains  neither  a  parameter  which  is  a  pair  nor  a  derivative  of  a 
pair,  The  computer  judges  whether  there  still  are  parameters 
■which  denote  -pairs  or  derivatives  of  pairs  in  location  26.0-30.15 
according  %o  register  0.10.  If  0.10  is  an  occupied  cell,  such 
parameters  ire  still  in  the  expression  and  the  repetition  of  cycles 
is  not  yet  finished.  If  0.10  is  an  enrpty  cell,  then  the  work  is 
.considered  finished  and  the  computer  is  stopped. 

The  basic  instructions  in  this  part  of  the  program  are:  the 
instructions  for  selection  of  the  element  form  location  1.0-5*15# 
for  checking  it  for  a  parameter  which  denotes  a  pair  or  the 
derivative  of  a  pair  and  storing  the  pair  or  derivative  of  the  pair 
selected  from  the  appropriate  location  massifs  6.0-10,15  or  11.0- 
23.15  in  location  63.II-63.I5#  the  instruction  for  recording  the 
element  in  location  26.0-pO.15  as  well  as  instructions  for  analysis 
of  the  elements  of  location  63.ll-63.i5  for  +,  comparing  previous 
elements  for  a  sign  of  an  elementary  function,  x,  :,  Y?  and  comparing 

the  subsequent  elements  for  x,  :,  . 

As  a  result  of  the  work  of  this  part  of  the  program  for 
the  above -examined  example,  we  will  obtain  the  following: 


26-0 

(1  :  *in 

26.1 

x)  1  cost 

26.2 

x  +  ex 

26.3 

(  ~  «. 

26.4 

x)-( 

26.5 

-a  )' 

26.6 

log  e«arct£ 

26.7 

x  ♦  1: 

26.8 

(1  *  X 

26.9 

*2). 

26.10 

e*(~ 

26.11 

a  •*) 

It  is  clear  that  by  storing  the  obtained  derivative  in  the 
cell  location  1.0-5.15  and  changing  the  same  program  again,  the 
second,  third,  etc.,  derivative  may  be  obtained. 

In  order  to  find  Taylor  coefficients,  to  solve  problems  in 
the  analysis  of  curves  and  surfaces,  to  find  maxima,  inflection 
points,  curvature  and  torsion,  etc.,  programs  are  needed  which,  in 
their  structure,  would  be  closely  related  to  the  program  analyzed 
here,  since  it  is  necessary  to  find  the  values  of  a  function  and  its 
derivative  at  a  given  point. 
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If  in  expanding  an.  expression  to  pairs  one  successively  finds  » 
the  numerical  Values  of  the  pairs,  then  the  numerical  value  of  the 
last  jpair  trill  be  the  numerical  value  of  the  whole  expression  at 
some  giveh  parameters. 

Thus,  the  first  part  of  the  program  for  automatic  differen¬ 
tiation  can  change  into  a  program  for  computation  of  the  value  of  a 
mathematical  expression  at  a  certain  point,  and  the  whole  program 
may  he  utilized  for  the  computation  of  derivatives  of  a  mathematical 
expression  at  that  point.  For  example,  in  order  to  obtain  the 
•numerical  value  of  k  derivative,  it  is  necessary  that  the  program 
he  repeated  k  times,  and  the  first  part  of  the  program  k  +  1  times. 

The  method  of  decomposing  an  expression  into  pairs  will  help 
to  solve  the  following  interesting  general  problem:  to  compose  a 
program  Op  which,  operating  with  the  commands  of  any  program  Pf 

that  computes  a  function  f (x),  would  develop  a  program  Ppf  that 

computes  a  function  F(f(x)),  where  Fg  is  a  certain  ope rater.  In 
particular,  by  F  one  can  mean  differentiation  or  integration. 

The  utilization  in  a  somewhat  different  manner  of  the  parts 
of  the  above-described  program  for  differentiation  of  mathematical  . 
expressions  will  allow  one  to  solve  the  problem  of  differentiating 
some  class  of  programs. 

In  actuality,  the  sequence  of  the  pairs  of  location  6.0-10.15 
is  nothing  else  but -a  program  for  computing  the  function  f  which 
was  originally  expressed  in  the  customary  mathematical  symbolism. 

The  two  last  parts  of  the  program  for  differentiation  allow  one  to 
obtain  a  mathematical  expression  for  the  derivative  of  the 
function  f.  It  is  evident  that,  if  the  obtained  expression  of  the 
derivative  is  expanded  into  elementary  pairs,  we  will  obtain  a 
program  for  computation  of  the  derivative  f For  this,  we  have 
only  to  apply  the  first  part  of  the  program  for  differentiation. 

Thus,  in  this  process  we  proceeded  fora  a  certain  program 
P we  applied  the  operator  0,-q  to  it  and  we  obtained  a  new  program 

Bof  —  a  program  for  the  computation  of  the  derivative  function 
of  f. 

As  a  conclusion,  we  will  cite  a  number  o^  examples  of  func¬ 
tions  and  their  derivatives  obtained  on  the  BESM  computer. 

Example  1 . 

f(x)  -  (arctg  In  (#“-  1)  +  ex)  t  (cos  Xs  -  *ln (**+  arc sin  *)) 
ftx)  -  (((1:  (l+(ln(*n-  l))8))  •  (l:{xB-l))‘*n-n.*x  +  sx*lne)* (eesfx8) - 
-sin(je*  ♦  aroiinx))  -((-sinfx8)* 

*x8  •  2  :x)  -  oob (x8*  areainx)  • 

• (x8*  3  :  x  *  I  :(1  -  x*)1^8  )). 

•(aretglnfx®-!)  ♦  e*)}  rfcoafx8)  - 
-  sln(x8+  arc  ■inf*)8 
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Example  2 

f(x)  -  (arctg  1ft  tgeos(x  +  arc  sin  *)) 

f'(x)  -  (1  :  (l+(lntg«>8(x*arcsinx))s)); 

♦(1  :  tgcos(x  ♦  arcsinx))*!  :(coscos  (x+  aresinx))®’ 

«(-sin(x  ♦  arc  sin*)  •  (1  +  1;  (1 
Example  3 

f(x)  -  (ln(x  *  lu(x  +  ln(x  +  la(  x  +ln*))))) 
f*(x)  “(1  :  (*  *  ln(x  +  1»(*  ♦  ln(x+  lw)))))* 

•  (1  +(1„*  (x  ♦  ln(x  +  ln(x+  lnx  ))))• 

• (l+(l: (x+ln(x+lnx)))* (l+(l:(x+lnx))*(l+l:x)))} . 


In  a  careful  analysis  of  the  derivative  formulas  obtained 
on  the  computer,  it  will  be  seen  that  further  simplifications  in 
them  are  possible.  Tills  problem  may  also  be  solved  by  means  of 
the  computer. 


Table  1 
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9 

IS 

h 
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1.6 

i 
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1.7 
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1.8 
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1 

1.9 
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1.11 

0 
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P 
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1.14 
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8 
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60 

1 

2 

3 

4 
6 
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sin  x 
la  6ft) 

«*x 

0  -  6,2 
#*6*3 
arctg  x 
6.4*6.6 


6 

9 

10 

11 
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4 


* 


li.O 

1 

z 

$ 

4 

5 

6 
? 

8 

9 

10 

11 

12 

13 
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(1 : 6.0 
)  t  ii.o 
& 

-11,2 
e*6.3 
11*3  4  log* 
1  :  (1 
+XY2 
) 

'  * 

11*4  *6.5  ♦ 

11*5  *6.4 
11*1  +  ii Je 
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Instruction  Diagram  1  of  the  program  for 
representation  of  a  mathematical  expres¬ 
sion  by  means  of  a  sequence  of  pairs. 


ftie  notation  »(*. y)  should  be  understood  as  follows: 
a  is  the  number  of  the  instruction,  £  is  the  number  of  the  instruc¬ 
tion  to  which  we  pass  in  case  of  positive  results  of  checking,  and 
Y  is  the  number  of  the  instruction  to  which  we  pass  in  the 
opposite  case.  Ihe  notation  <*U  -  t)  .means  that  after  the  instruc¬ 
tion  I'Iq,  one  should  pass  to  the  instruction  Up,  and  then  to  -the 

instruction  N  y  . 


1(2,2) 

2(1,5) 

3(1,4) 


* 


4(1,5) 
5(28,6) 
6(16,7) 
7(8, io; 
8(9,18 
9(1,20, 
10(11,13) 
11(12,17) 


The  selection  of  the  next  symbol 
Checking  for  the  sign  of  an  open  bracket 
Checking  for  the  sign  of  an  elementary  function  and 
parameter 

Checking  for  the  sign  of  a  power 

Checking  for  a  multiplication  and  division  sign 

Checking  for  an  addition  sign 

Checking  for  a  subtraction  sign 

Checking  the  preceeding  sign  for  an  open  bracket 

Checking  the  subsequent  sign  for  an  open  bracket. 

Checking  for  the  sign  of  a  closed  bracket 

Checking  upwards  over  one  for  the  sign  of  an  open  bracket. 
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12(1.1) 

13(14,15) 

14(0,0) 

15(0,0) 

16(1,21) 

17(11,11) 

18(9,16) 

19(1,26) 

► 

20(1,1) 

21(22,23) 

22(16,16) 

23(24,22) 

24(25,19) 

25(2,2) 

26(28,28) 

27(24,26) 

28(1,27) 


1(2,2) 

2(3,4) 

3(0,0) 

4(5,29) 

5(6,6) 

6(7,11) 

7(8-9) 

8 

9(10,10) 

10(1,1) 


11(7,12) 

12(7,13) 

13(14,16) 


Destruction  of  brackets 

Checking  for  the  sign  of  the  end  of  a  formula  (empty) 

Stop  (correct) 

Stop  (incorrect) 

Checking  over  one  symbol  upwards  for  the  sign  of  an 
open  bracket 

Recording  a  pair  in  location  6.x.  Exchange  of  variables 
Checking  the  proceeding  sing  for  the  sign  of  a  power 
Checking  for  multiplication  and  division  the  last  sign 
selected  from  the  expression 

Recording  the  pair  in  location  6.x.  Exchange  of  variables 
Checking  over  one  symbol  upwards  for  a  sign  of  an 
elementary  function 

Recording  the  pair  in  location  6.x.  Exchange  of  variables 
Checking  for  the  sign  of  subtraction 
Checking  the  sign  before  a  minus  for  a  power 
Recording  a  pair  in  the  location  6.x.  Exchange  of 
variables 

Recording  a  pair  in  the  location  6.x.  Exchange  of 
variables 

To  check  over  one  symbol  upwards  for  a  subtraction 
sign 

Checking  over  one  symbol  upwards  for  signs  of  addition 
and  open  bracket. 

Instruction  Diagram  2  of  the  program  for  obtaining 
derivatives  of  elementary  pairs 

Selection  of  the  pair  from  location  6.0-10.15 
Checking  for  emptiness  in  the  cells  of  location 
6.0-10.15 
Print  shop 

Checking  for  pairs  which  consist  of  two  symbols 
Formation  of  derivatives  from  the  parameters  which 
enter  the  given  pair 

Checking  the  first  symbol  of  the  pair  for  sin 
Formation  of  a  derivative .from  sin.  v,  cos  v,  exp  v. 
Recording  derivatives  in  general  form  into  the 
location  of  cells  59*0-59.4 

Simplification  of  the  expression  of  derivatives  from 
sin  v,  cos  v,  exp  v 

Rs cording  the  final  expression  of  the  derivative  of  a 
pair  in  location  11.0-25.15  and  recording  the  ordinal 
number  of  the  cell  in  which  the  derivative  originates 
in  the  last  8  columns  of  the  cells  of  location 
6. 0-10. 15  which  correspond  to  the  given  pair 

Checking  the  1st  symbol  of  the  pair  for  cos 

Checking  the  1st  symbol  of  the  pair  for  exp. 

Checking  the  1st  symbol  of  the  pair  for  In 


*14(8-15) 

15(10,10) 

16(17,19) 

17(8-18) 

l8(io,io) 

19(20,21) 

20(8-18) 

*21(22,24) 

22(8-25) 

£5(10,10) 

24(22,25) 

25(26,28) 

26(8-27) 

27(10,10) 

28(0,0) 

29(30,50) 

30(31,33) 

31(8-32) 

32(10,10) 

33(34,36) 

34(8-35) 

35(10,10) 

36(37,39) 
37 (8-38) 
38(10,10) 

39(8-40) 

40(10,10) 


1(2,2) 

2(3,3) 

3(7,4) 

4(5,22) 

5(12,12) 


Formation  of  a  derivative  of  In  v 
Simplification  of  the  expression  for  a  derivative 
of  In  v 

Checking  the  first  symbol  of  the  pair  for  tg 
Formation  of  the  derivative  of  tg  v 
Simplification  of  the  expression  of  the  derivative 
of  tg  V  and  ctg  v 

Checking  the  first  symbol  of  the  pair  for  ctg  v 
Formation  of  a  derivative  from  ctg  v 
Checking  of  the  first  symbol  for  arc  sin 
Formation  of  a  derivative  of  arc  sin  v  and  arc 
cos  v 

Simplification  of  the  expression  for  a  derivative  of 
arc  sin  v  and  arc  cos  v 

Checking  the  first  symbol  of  the  pair  for  arc  cos 
Checking  the  first  symbol  of  the  pair  for  arctg 
Formation  of  a  derivative  of  arctg  v 
Simplification  of  the  expression  of  the  derivative 
of  arctg  v 
Stop 

Formation  of  derivatives  from  the  parameters  which 

enter  into  the  given  pair 

Checking  for  multiplication  sign 

Obtaining  a  derivative  from  u*v 

Simplification  of  the  expression  for  a  derivative 

of  the  product 

Checking  for  a  division  sign 

Formation  of  a  derivative  from  u:v 

Sircplifi cation  of  the  expression  for  the  derivative 

of  the  quotient 

Checking  for  the  sign  of  a  power 
Formation  of  a  derivative  from  uv 

Simplification  of  the  expression  for  a  derivative  of 
the  power 

Formation  of  a  derivative  of  u±v 

Simplification  of  the  expression  of  the  derivative  of 
the  sum  and  difference 

Instruction  Diagram  3  of  the  program  for 
synthesis  of  the  derivative  formula 

To  select  the  derivative  of  the  last  pair  and  record 
it  in  location  1.0-5. 15 

To  select  the  next  element  from  the  location  1.0-5.15 
To  check  it  for  being  empty 

To  check  it  for  "u"  or  "u'",  i.e.,  is  the  element  a 
pair  or  the  derivative  of  a  pair? 

To  select  the  value  u  (from  location  6.0-10.15)  or 
u'  (from  location  11.0-25.15)  and  record  it  in 


5(0,0) 

7(6,21) 

3(16,18) 

9(16,17) 

1-0(9,15) 

11(16,10) 

12(17,11) 

15(19,8) 

14(16,16) 

15(16,14) 

16(2,2) 

17(2,2) 

18(16,17) 

19(16,20) 

20(16,17) 

21(2,2) 

22(2,2) 


location  65.11-63*15 

The  printing  of  results  (location  6.0-10. 15,  11.0-25.15, 
26.0-30.15).  Stop 

To  check  for  the  end  of  substitution,  i.e.,  for  the 
absence  in  26.0-30.15  of  u  and  u' 

To  check  the  element  before  the  selected  element  in 
location  I.O-5.I5  for  a  sign  of  an  elementary  function 
To  check  the  element  next  after  the  selected  element 
in  location  1.0-5.15  for  a  sign  of  a  power 
To  check  location  63.U-63.i5  for  two  elements 
To  check  the  first  element  in  63. 11  for  a  subtraction 
sign 

To  check  the  location  63. 11-63*15  for  one  element 
To  check  the  location  63.ll-63.i5  for  signs  of  addition 
and  subtraction 

To  check  the  element  following  after  the  selected 
element  in  location.  1.0-5.15  for  the  sign  of  a  power 
To  check  the  element  next  after  the  selected  element 
in  location  1.0-5.15  for  a  division  sign 
To  record  the  content  of  cells  63.H-63.i5  in  location 
26.0-30.15  with  brackets 

To  record  the  content  of  cells  63. 11-63.15  in  location 
26.0-30. 15  without  brackets 

To  check  the  element  next  after  the  selected  element 
in  location  I.O-5.I5  for  a  sign  of  a  power 
To  check  the  element  before  the  selected  element  in 
location  1.0-5.15  for  a.  sign  of  multipH cation,  division 
degree  or  elementary  function 

To  check  the  element  following  after  the  selected 
element  in  location  1.0-5*15  for  signs  of  multiplica¬ 
tion  or  division 

To  transmit  the  whole  phrase  from  location  26.0-30.15 
into  location  1.0-5.15*  To  clear  location  26.0-30. 15 
To  record  the  element  selected  from  location 
1.0-5.15  into  location  26.0-30.15  without  a  bracket 
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